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Abstract 

Ice thickness is an important variable for climate scientists and 
is still an unsolved problem for satellite remote sensing specialists. 
There has been some success detecting the thickness of thin ice from 
microwave radiometers, and with this in mind this study attempts to 
model the thickness-radiance relation of sea ice at frequencies employed 
by the Soil Moisture and Ocean Sahnity (SMOS) radiometer and the 
Advanced Microwave Scanning Radiometer (AMSR): between 1.4 and 
89 GHz. In the first part of the study, the salinity of the ice was deter- 
mined by a pair of empirical relationships, while the temperature was 
determined by a thermodynamic model. Because the thermodynamic 
model can be used as a simple ice growth model, in this, second part, 
the salinities are determined by the growth model. Because the model 
uses two, constant-weather scenarios representing two extremes ("fall 
freeze- up" and "winter cold snap"), brine expulsion is modelled with 
a single correction-step founded on mass conservation. The growth 
model generates realistic salinity profiles, however it over-estimates the 
bulk salinity because gravity drainage is not accounted for. The results 
suggest that the formation of "skim" on the ice surface is important 
in determining the radiance signature of thin ice, especially at lower 
frequencies, while scattering is important mainly at higher frequencies 
but at all ice thicknesses. 



1 Introduction 



This study is a follow-up to a similar study ( Mills and Hevgster . 201 ll ) in 
which the relationship between sea ice brightness temperature at microwave 
frequencies and ice thickness was examined. As noted in the previous study, 
several other works have demonstrated a relationship between ice thick- 



ness and SSM/I o r AMSR-E brightness temperatures (IMartin et al.l . l2004l : 



Naoki et al.l . l2008l : iKwok et al.l . 120071 : iHwang et al.l . 120071 ) with at least one 



attempt to use this relationship for the purpose of ice thickness retrieval 
(jMartin et al.l . 12004 ) . To understand this relationship, determine it's most 
likely functional form or forms and to set definite bounds on it would be of 
great value to ice remote sensing efforts. In particular, satellite retrieval of 
sea ice thickness is still an unsolved problemn, yet is invaluable for climate 
scientists as it can help determine a number of important flux variables in 
the arctic regions such as heat and salt and fresh water. 

In the first part, parameterised empirical relationships were used to relate 
bulk salinity to ice thickness and to determine the vertical salinity profile 
within the ice. In this second part, we will use a thermodynamic ice growth 
model, based primarily on first principles, to determine the salinity of the 
ice. 

In Section [2] we will describe in detail the thermodynamic ice growth 
model, while in Section [3] we will present the results of the that model. 
Finally, in Section H] we will apply the modelled profiles to the microwave 
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Table 1: Thermodynamic model parameters for the two weather scenarios 





Scenario 


Parameter 


Fall freeze-up Winter cold snap 


Wind speed [m/s] 
Air temperature [K] 
Relative humidity 
Cloud-cover 
Insolation [W/m^] 


2. 10. 
270. 260. 
0.5 0.1 
0.5 0.1 
50. 0. 



emissivity simulation. All other parameters, in particular the emissivity 
model set up, will be the same as in the previous study and so will not be 
described. 



2 Ice growth model 



In the first part of the study ( Mills and Hevgster . 201ll ) a simple ice growth 



model based on constant weather scenarios was described. Here we apply 
the same model but with an extra step modelling brine expulsion added. 
These brine expulsion processes result from cooling of the ice sheet in the 
upper layers: because ice has a lower density than brine, as the ice grows 
downward, brine in the uppe r layers will f r eeze generating pressure in the 



remaining, more saline brine (jTucker et al.l . Il992l ). 



Since we are using two, constant-weather scenarios-see Table [T]-representin^ 
two possible extremes, there is no need to created an iterative, time-based 
prognostic scheme, rather we take only ice thickness, h, in favour of time, t, 
as the dependent variable in the bulk of the model. The salinity profile is 
first modelled assuming no brine expulsion, then there is a "brine correction" 
step, in which mass conservation determines how much brine is lost based 
on the difference in temperature between the ice layer at first formation and 
at the final ice thickness. 

The following equation relates the ice surface temperature to net heat 
flux: 

hQ* = k{T„ - Ts) (1) 

where r„ is the water temperature which is assumed to be constant at 
freezing (approximately -1.9° C at a water salinity of 35 psu), Tg is surface 
temperature and k is the thermal conductivity of the ice. The net heat flux 
comprises the following components, with functional dependencies supplied: 

Q* = QE[e{Ts)] + Qh(Ts) + Qsw{T,^) + Qlw (2) 
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The terms on the RHS are, from left to right: latent heat, sensible heat, 
longwave and shortwave; e(T) is the saturation vapour pressure. The first 
two terms are approximated with simple parameterisations while the long- 
wave flux is based on the Stefan-Boltzmann law. The shortwave flux is 
calculated primarily from geometric considerations based on the position of 
the Earth relative to the Sun. The following inputs are required for the 
model: surface- wind speed, humidity, air temperatu re and density (o r pres- 
sure), cloud cover and date and time or insolation. ( Cox and Weekj . Il988l : 
Drucker et all . l2003l : IYu and Lindsavl . l2003l ) Inputs can be supplied to give a 
picture of the general weather conditions. For instance, fall freeze-up might 
be characterized by relatively mild temperatures, low winds, high humid- 
ity, high cloud cover and moderate insolation. By contrast, a winter cold 
snap would be characterized by low temperature, high winds, low humidity, 
clear conditions and little to no insolation. See Table [TJ Equations ([1]) and 
( 21) are solved with a numerical root-finding algorithm, specifically bisection 



Press et al.l . Il992l 'l. 



The model can also be used as a crude ice growth model. The rate of 
ice growth is: 

where L is the latent heat of fusion for water and pi ~ 0.917 — 1.403 x 10 
is the density of pure ice, which is a function of temperature. Empirical 
equations f or determining the init ial br ine entrapment in sea ice have been 
derived by ICox and Weeksl (|l988l ^ and Nakawo and Sinhal (|l98ll ) and take 
the form: 

S = Sofig) (4) 

where Sq is the salinity of the parent water and / is an empirical function 
of ice growth rate. 

As shown in the previous report, if the salinities are left untouched, 
even as they move further from the bottom of the ice sheet, the model 
considerably over-estimates the bulk salinities. Using conservation of mass, 
we can correct for this in each layer based on the temperature differences 
between the thinner ice, when the layer is still forming, and the thicker ice as 
more layers have been added to the bottom. In the following, we go through 
all the relevant relationships needed to form the brine expulsion-correction 
model. 

While we will not need it, the time, t, can be written as a function of ice 
thickness: 
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^\h')dh' 



(5) 



since the growth rate, g, is also a function of ice thickness. Because we 
assume constant weather conditions, temperature for a given point in the 
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ice sheet is a monotonic function of ice age or thickness, hence there is no 
reason to construct a dynamical time series. 

Due to freezing-point depression, the sahnity of brine inclusions, Sb, is a 
function of temperature only. That is, greater salinity in the brine will lower 
it's freezing point. As the brine cools, the water in it will freeze, increasing 
its salinity. Thus it is maintained at its free zing point at all tim es. We use 
the linear, piece- wise continuous model from Ulabv et al. ( 19861 ). 



The total density at a given point in the ice is: 

P = Vi,Pb + (1 - (6) 

where Vh is the r e lative brine volume, pb ~ 1 + 0.00085'b is the brine density 
(|Cox and Weeksl . Il983l 'l. Note that this neglects any air pockets that are 



usually included in the ice in addition to brine. These are particularly 
prevalent in older ice, but omitting them considerably simplifies the analysis. 
The ice salinity is further related to the brine salinity as follows: 

S = 5b^ (7) 
From equations ([6]) and ([7]) we can solve for the relative brine volume, V\y: 

= ^ (8) 

We assume that the ice is in thermodynamic equilibrium and the thermal 
conductivity is roughly constant throughout. Thus, the temperature of the 
ice sheet at depth, z, and thickness, /i, is given by interpolating between the 
surface temperature, Tg and the water temperature, T^: 

T{h, z) = T„ + [T,{h)-T^]z/h (9) 

From this we can calculate the density of the ice sheet as a function of ice 
thickness and ice depth: p{h, z). The relative change in volume, 5v, is given 
as: 

Sv = {p/po - 1) (10) 

where po{h) = p{h, h) is the density of the ice layer as it is first forming. 
Since the amount of mass must remain constant, and since we assume there 
are no cavities in the ice, so must the volume remain constant, therefore we 
can determine how much material is ejected from the ice sheet. We assume 
that only brine is ejected, provided there is enough, thus: 

5m = h^""'^ - (11) 

\phVh + Pi{5v -Vi,); 5v>Vi, 
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2 4 6 610 b 10 15 20 

S [psu] S [psu] 

Figure 1: Modelled salinity profiles for different ice thicknesses and two 
weather scenarios: "fall freeze-up" and "winter cold-snap"; see Table [H 



is the relative change in mass or relative mass loss. It follows that the 
corrected brine volume is given as: 



[0; Sv>V^ 
from which we can solve for the corrected salinity using Equation ([8|). 

3 Growth model results 

Results for the "fall freeze-up" and "winter cold-snap" scenarios are shown 
in Figures [U^a) and (b), respectively. Qualitatively, at least, these salinity 
profiles are quite realistic: thin-ice profiles for the c old-wea t her s cenario 



show the characteristic, 'C'-shape typical of young ice (lEickenl. Il992l 'l. while 
the thicker, older ice has a flatter profile ( Tucker et al. . 19921 ). For the 



warm- weather sce nario, the profiles are altogether fiatter, which is typical 
of low-salinity ice (|Granskog et al.l . l2006l ). 



Bulk salinity as a function of thickness are shown in Figure [21 while 
thickness as a function of time is shown in Figure [3l There are two versions 
of the growth model salinities: "with skim" and "without skim", both of 
which deserve some explanation. The brine that has been ejected has to go 
somewhere. In the first case, we assume that exactly half of it collects on 
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Salinity — thickness relation 




Figure 2: Modelled versus empirical salinity-thicknes s relat ionships. Solid 
shading shows the residuals for the Cox and Weeks ( 19741 ) models. The 
areas in between the two weather scenarios are hatched. See text for an 
explanation of the difference between "with skim" and "without skim." 
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t [days] 

Figure 3: Ice thickness as a function of time for the two contant- weather 
scenarios. 



top of the ice sheet as "skim" , while in the second case, we assume that all 
of it is ejected downward, into the ocean. 

We need to calculate both the salinity of the skim, and its thickness for 
further use in the emissivity model. Following Equation d?]), the lost salinity 
in each layer is: 

5,ost = |^''''''/'' (13) 
[ShPhVh/p; 6v>Vb 

The total mass loss is given as: 

Srritotai = / Smdz (14) 



thus the salinity of the skim is calculated: 



1 



h 



Sskim = T / Siost^mdz (15) 

Assuming that the skim is at the temperature of the top layer, we can 
compute the density of the skim, Pskinn from Equations ([S]), then ([U]), from 
which follows its thickness: 

C>2;skim = ^3 (16) 

■^Pskim 
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Figure 4: Thickness and salinity of "skim" layer. 



Properties of the modelled skim layer are shown in Figured 

Note that even if all the ejected brine is lost from the ice sheet, the 
model still over-predicts the salinities in relation to the Cox and Weeks 
model, especially at deeper ice thicknesses. The reason for the discrepency 
is that brine drainage has not been included in the growth model: as the ice 
ages, brine will drain through vertical channels. Still, we hope that feeding 
these results to the emissivity models will nonetheless provide some insight. 
For one thing, the salinity profile is different; also, there are more emissivity 
variations in thinner ice, where the model is more accurate. 



4 Results 

Here the final brightness temperature-thickness relations are presented. In 
this second part, the error bounds h ave also been calc ulated using a Monte 
Carlo or "boot-strapping" method ( Press et al. . 19921 ) by first varying the 



constant weather scenario between the two extremes presented in Table [T] 
to produce 100 sets of profiles for different ice thicknesses, then using the 
error bounds in the effective permittivity model to vary the effective per- 
mittivities over 100 additional trials, generating 10000 trials total. Uniform 
deviates ( "white" noise) are used to vary the weather parameters while nor- 
mal deviates (Gaussian noise) are used to vary the permittivities. 

In the first set of graphs, Figure [5l neither the skim layer nor the scatter- 
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(a) Brightness temperature vs. ice thickness 
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Figure 5: Modelled brightness temperature as a function of ice thickness 
based on modelled ice profiles, not including skim. Solid lines show model 
results neglecting scattering, while dotted lines show the computed scatter- 
ing component. 
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(a) Brightness temperature vs. ice thickness 
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'b) Polarization differenoe vs. ice thickness 
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Figure 6: Modelled brightness temperature as a function of ice thickness 
based on modelled ice profiles, including skim. Solid lines show model re- 
sults neglecting scattering, while dotted lines show the computed scattering 
component. 
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[a) Brightness temperature vs. ice tiiickness 
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Figure 7: Modelled brightness temperature as a function of ice thickness 
based on modelled ice profiles, not including skim. Scattering has been 
included. 
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[a) Brightness temperature vs. ice thickness 
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(b) Polarization difference vs. ice thicl<ness 
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Figure 8: Modelled brightness temperature as a function of ice thickness 
based on modelled ice profiles, including skim. Scattering has been included. 
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ing component has been included in the emissivity calculations. The scatter- 
ing components, however, have been graphed alongside the non-scattering 
components. The second set of graphs, Figure [U is the same as the first, 
except the skim layer has been included. In the third set of graphs, Figure 
[71 the scattering and non-scattering components have been added together 
for the simulations not including skim. The fourth set of graphs. Figure [8l 
is the same as the third, except the skim layer has been included. 

One of the most striking things about these results is how uniform the 
brightness temperatures are in are in Figure O both for different ice thick- 
nesses and for different frequencies. Only the lowest frequency shows a 
slight dip at the thinnest ice thicknesses because of transparency: the water 
"shows through." Even with the addition of skim in FigureEl the brightness 
temperatures are still surprisingly uniform. Skim is one factor that is rarely 
included ice growth or emissivity models. Because of its high liquid and 
high saline content, it will be highly reflective. Thus its presence will tend 
to lower the brightness temperature of the ice because it reflects the much 
cooler (as little as 3 K on clear days) sky, as well as increase the polariza- 
tion difference because of its higher real permittivity. These results show 
a fairly uniform depression/increase in brightness temper ature/polarization 
difference, but because the modelled skim does not drain as it would on real 
ice, the actual tendency would apply only to new and thin ice. 

Because the both real imaginary permittivity increa se more with brine 
volume the lower the frequency (see Figure 10 from Part I (jMills and Hevgsteii . 



201 ih ). the effect of the skim is more pronounced at lower frequencies. Thus, 



in lower frequencies, the characteristic lowered brightness temperature and 
raised polarization difference of new and thin ice is caused by the accumula- 
tion of skim. In higher frequencies, this is caused by scattering. The results 
here show that scattering increases the variability of the signal especially at 
higher frequencies and especially at higher ice thicknesses. Note that skim 
also increases the variability, but more so at lower frequencies. In real ice we 
would expect the effect of scattering to be more pronounced in thin rather 
than thick ice, and this does show up somewhat in the model results, but 
only when skim is included. 



5 Summary and conclusions 

A simple thermodynamic model was used to simulate ice growth for con- 
stant weather scenarios ranging between two extremes. This model takes 
into account both heat conduction and brine expulsion without the need 
to construct a dynamical time series, but neglects brine drainage. The 
resulting temperature and salinity profiles are fed to the Microwave Emis- 
sion Model of Layered Snowpack (MEMLS) to determine the functional 
dependence of microwave radiance on ice thickness. The resulting salinity 
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profiles and consequent salinity-thickness relation, as well as the simulated 
brightness-temperature-thickness relations, are quite realistic. The follow- 
ing observations may be noted: 1. Without a layer of "skim" collected on 
top of the ice and neglecting scattering, the brightness temperatures of the 
ice are surprisingly uniform, both over different frequencies and different 
ice thicknesses. Thus, variations in emissivity must be attributed to other, 
secondary factors which we describe below. 

The presence of skim, a semi-liquid collection of brine and ice that accu- 
mulates on top of ice sheets due to brine expulsion, is neglected in most ice 
growth and emissivity models. Here the skim has been modelled as simply 
the aggregate of exactly half the brine that has been expelled from the ice 
sheet. 2. Because of its high opacity and high real permittivity, the pres- 
ence of skim will tend to lower the brightness temperature and increase the 
polarization difference. This effect is more pronounced at lower frequencies 
and at thinner ice thickness, but the latter effect is not present because the 
growth model neglects brine drainage. 3. Thus, the lower brightness tem- 
perature and higher polarization difference characteristic of new and thin 
ice is primarily caused by the presence of skim. 

4. Finally, scattering will have a strong effect at higher frequencies. 5. 
Skim in particular will scatter strongly. 6. Thus, lower brightness tem- 
perature and higher polarization in thin and new ice and lower brightness 
temperature and higher polarization difference in general at higher frequen- 
cies is caused by scattering, especially by the skim. 7. Both the presence of 
skim and scattering will tend to increase the variability of the signal. 

Both ice growth and the interaction of sea ice with electromagnatic radi- 
ation are extremely complex processes with many different facets. They are 
very difficult to model effectively. Moreover, since all of: the microstructural 
properties, the macroscopic configuration, the relative composition, will af- 
fect the emitted microwave signal, and all these can change with the pre- 
vailing weather conditions, ice growth and ice emissivity cannot be treated 



in isolation. Figures 1 and 2 in Part I (jMills and Hevgsterj . |2011| ) give some 
idea of the complexity of the problem. 

We believe the contents of this and the previous report represent an 
advance in the state-of-the-art, however there are still many limitations in 
the current effort. 1. The ice growth model over-estimates salinity and 
the size of the skim layer at high ice thicknesses because it neglects brine- 
drainage. 2. Scattering is too high and too variable at high ice thicknesses. 
The scattering model is fairly crude and likely innaccurate, not least of which 
because the correlation length are estimated by a simple, ad-hoc assumption. 
3. The complex permittivities do not take into account ice microstructural 
properties, they do not have a comprehensive physical explanation and they 
are not validated at higher frequencies. As a first step, there are many new 
directions that the semi-empirical permittivity model could be taken in. A 
Debye relaxation curve for different brine volumes, for instance, is suggested 
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by the strong, simple dependence on frequency. 
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